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The improvement of the techniques for cooling atoms has been boosting advances 
in physics for several years [D El El i. A remarkable example is the hnding of 
experimental evidences of the so-called Efimov states, that were predicted by V. Efimov 
in the seventies [5]. Using few-body techniques, he showed that the energy levels of 
three identical bosons accumulate with a geometrical separation between them when 
the binding energy of the pairs goes to zero. More than thirty years passed until a 
signature of such states was found in loss experiments in 2006 |6]. Although being a 
true breakthrough, this was still an indirect measurement and another nine years were 
necessary until this effect was found in a direct manner in experiments with '^He |7]. 

This impressive technological advance also opened new avenues as the quasi-2D 
realization of single-species systems [HI El EQ], two-component gases mills], mixed- 
species systems [131 113] and even heteronuclear diatomic molecules [I5] in ultracold 
atomic traps. On the theoretical side, it has been shown that three-identical weakly 
attractive bosons in 2D have only one two-body and two three-body bound states [16] 
and further studies have also shown that the Efimov effect does not happen when 
the system is restricted to lower dimensions mm- Recent studies have focused on 
mixed-species two-dimensional systems [201 EHl [21] ; which, despite of not supporting the 
Efimov states, can have a rich energy spectrum in the case of highly mass-imbalanced 
systems [22l [23] . 

While most of these calculations were done in vacuum or disregarding the medium, 
bound states can be drastically affected depending on the surrounding medium in which 
the system is immersed. If one thinks about a single impurity immersed in a bath of 
different atoms, its motion distorts the medium, which in turn acts back on the impurity. 
The final state depends on the interaction strength between the different species, where 
for weak enough attraction the impurity behaves as a free particle and for strong enough 
attraction it forms a bound state with a particle from the sea. On the other hand, for 
intermediate attraction the impurity is dressed by the particles from the sea and becomes 
a polaron, namely, a quasi-particle whose basic properties as effective mass and charge 
differs from the ones of the original impurity. These different properties result from the 
dressing of the impurity by excitations of the background medium. The original idea of 
polaron dates back to Landau [23] and a discussion about it in the context of ultracold 
atomic gases can be found in [25] . 

It is possible to simulate polaronic systems with cold atomic gases by mixing a few 
atoms of a specihc type to a bath of atoms of a different kind [26] (different masses 
or spin states, for example) even when the system is restricted to two dimensions 
and the majority atoms are fermions, in which case the quasi-particle is known as a 
Fermi polaron. Fermi gases have been successfully trapped to quasi-2D geometries 
and the Fermi polaron was already observed through Radio-Frequency Spectroscopy 
techniques [27] EH] EH] EO]. The fate of a single impurity immersed in a two- 
dimensional Fermi sea has been studied in recent years both experimentally [3l] and 
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theoretically [521 ESI EH ESI ES], where the importance of the particle-hole fluctuations 
in 2D systems was emphasized. 

In real systems it can be very hard to have only a single particle propagating in the 
medium and it is therefore important to understand how the presence of extra impurities 
affect the system. In this work we take a hrst step in this direction and consider the case 
of two identical impurities immersed in a Fermi sea. The study applies to a broad range 
of masses and both for interacting and non-interacting impurities. In both cases, the two 
impurities can attach to a fermion from the sea and form three-body bound states. We 
show how the energy of the states are shifted depending on the Fermi momentum kp and 
also how these states decay into the two- and three-body continuum. The corrections 
due to the fluctuations in the sea are considered when the impurities are not allowed to 
interact with each other. We show that in spite of the strong influence both the two- 
and three-body systems, the corrections from fluctuations in the Fermi sea have small 
influence on the point where the three-body bound state cease to exist. 

This paper is organized as follows: the influence of the background environment 
in the state of an impurity bound to a fermion from the sea is discussed in Sec. [2] 
and the derivation of the equations which consider corrections to the two-body system 
due to the particle-hole fluctuations in the sea are presented in Sec. |3l Sec. 0] brings 
a discussion about the dependence of the two- and three-body thresholds in the Fermi 
momentum. Three-body systems composed of two identical interacting impurities and a 
fermion from the sea are studied in Sec. [5l Corrections to the three-body states induced 
by the fluctuations in the sea are discussed in Sec. [6] for two non-interacting particles, 
where it is also argued how the corrections would affect the three-body system when 
the impurities do interact. Discussion and outlook are in Sec. [7] and technical details 
are given in 


Appendix A and Appendix B 


2. Bound and virtual states of the two-body system 


We consider an atom labeled a of mass nia interacting with a fermion b of mass mb that 
is on top of a background Fermi sea of momentum kp. The interaction is described for 
the one term separable potential V = \ \x) (x|, it is attractive (A < 0) and has the form 
factor gfp) = (p] x) = 1 because we use a Dirac delta potential. The matrix elements 
of the two-body T-matrix [3711I9] under the influence of the Pauli blocking effects are 
calculated as 


t-\E^) = A-' 



p- 

^ rria ^ 




le 


( 1 ) 


where p is the relative momentum, q is the center-of-mass (CM) momentum of the pair 
with respect to the Fermi sea, the reduced mass is mat = mamb/{ma + mb) and E 2 is the 
internal (or binding) energy of the pair. The ©—function ensures that the fermion is not 
inside the Fermi sea. The divergence in ([1]) is treated with the subtraction method |38j . 
where the strength of the potential A is connected to the binding energy of the pair \Eab\ 
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through X ^ = — J d‘^p{\Eab\ + p^/2mah) Performing the integral in ([T]) yields 

/ 

= -2'Kmab'hn 


[a Ilf 

1/2 

[2E 

1/2 

2\Eab 



V 

- - E2 


2mab 2m2 


2\E. 


ab 


( 2 ) 


Notice that if the background Fermi sea is absent, kp = 0 and ([2]) recovers the two-body 
T-matrix elements of an ab system in vacuum, namely 





( 3 ) 


It is possible to identify three different thresholds in (|2]). As discussed in [3l] for 
= '<^b, two of them come from the branch cut of the square root and are given by 

E^ = ■ Thfi third one arises from the branch cut of the logarithm 

and reads = 0 for g > ^^kp. 

The two-body system only supports bound states when (i) there is a pole in the 
two-body T-matrix and {ii) the energy where the pole appears is below the lowest 
threshold, namely E 2 < Efj^ or E 2 < E^j^ for q > ff^kp. Setting ([2]) equal to zero, the 
two-body energies corresponding to bound states are described by the expression 


E 2 — —\Eab 


kp 

2mab 1 


1 


m-ab ’ 

2m| \Eab\ 


( 4 ) 


which tends to the energy of the pair in vacuum when the Fermi sea is not relevant, 
i.e., E 2 —)■ —\Eab\ for q/kp —)■ 00 . This expression agrees with the ones calculated using 
many-body techniques, variational methods etc. |32l |33l [SU |39l l36] . 

The internal two-body energy given in (jl]) always satishes the condition (zi), but 
surprisingly it does not satisfy the condition (i) for all possible combinations of momenta 
q and kp. There is a strong competition between binding, the Fermi sea and the CM 
momentum such that the pair is unbound for kp > SrOab \ Eab \ when the value of these 
parameters fulhll the relation 

(^F - \Jkp- 8mab\Eab\^ <q< (^kp + \J kf - Smab\Eab^ -(5) 

Although the analytical expression (jl]) is a solution of t~^{E 2 ) = 0, numerical 
calculations from ([T]) also do not indicate two-body bound states in the region given 
in ([5]), which has its boundaries framed by the vertical arrows in hgure 

The absence of bound states in the region dehned by ([5]) was already pointed out 
in [M] and can also be seen in [ID] for the spin-balanced system and the remaining 
question is to know what happens to the state in this region. It is shown in hgure [T] 
that for kp > ^/^mab\Eab\ the state touches the lowest threshold {Ef^) in two points. 


1(b) 
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strongly suggesting that it enters through the two-body cut in one point and conies 
back in another one. This cut is determined exclusively by the square-root inside 
the logarithm in ([2]) and the logarithm itself does not contribute to the cut since its 
argument is positive and non-zero throughout the regions where E 2 < and E 2 < E^j^ 
for q > ^kp. 


The analytic continuation of the matrix elements ([2]) is labeled and found 
by flipping the sign in front of the square-root inside the logarithm in ([2]), since the 
logarithm itself does not contribute to the cut. Next, this expression is used to calculate 
E 2 as a solution of t~^{E 2 ) = 0, which gives precisely (jl]) as the result. Notice that 
now the solution still satishes condition (ii), but only satishes condition {i) in the region 
given in ([5]) and not outside this region as before. These two-body states have entered 
through the cut and are in the second energy-sheet, therefore they are virtual states. 


These states are shown in hgure 1(b) 



(cl) kp SlTlab I tt/ab \ (h) kp > \/ ^niah \ t^ah \ 


Figure 1. Two-body internal energy, E 2 /\Eab\, as function of the CM momentum 
q/\/ma\Eab\- The analytic expression for the two-body energy in (jH) perfectly matches 
numerical calculation for both bound and virtual states. |l(a)| Only bound states are 
present for kp < \/8mab\Eab\- |l(b)[ Both bound and virtual states are present for 
kp > \/8mab\Eab\- The vertical arrows are the extremes of jS]). 

The matrix elements of the two-body T-matrix as given in ([2]) may lead to 
inconsistent results in the scattering region {E 2 > Efffj, where the imaginary part of the 
energy becomes important. For instance, setting E2 E2 ie in ([2]), expanding the 
terms inside the square root and collecting energies and e’s together makes it possible to 
see that the imaginary part becomes null in some cases. If one disregards for a moment 
the small imaginary term that should appear inside the square roots and considers that 
the imaginary contribution due to these terms come only from the arguments inside 
it, then three cases have to be considered, namely, E 2 < Efj^, Efy^ < E 2 < E'^^ and 
E2 > El. 

When the energy crosses from E 2 < E^j^ to E 2 > E^f^, the analytic extension of the 
square root has to be properly taken in account. In other words, its imaginary part has 
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to follow the sign of the ie term that appears outside the square roots in ([2]), meaning 
that [l/2mab {kp - rnahjma qf - -t -i [E 2 - l/2mab {kp - ruab/ma for 

E 2 > Efi^. The same argument is taken in the transition from Eff^ < E 2 < Ef^ to 
E 2 > and therefore, the matrix elements are written as 

' g{E2,kp,q) + ^^{kl 




{E 2 ) = -2pmab hi 


^lb„2 

T7,2 y 


Eo 


le 


2\E, 


ab 


■ ( 6 ) 


with g{E 2 , kp, q) given by 


for E 2 < E^^, 


-^(kp-^q] -E2 
2mab V 


E 2 - - - kp - q 

2mab V 'TUa 


1/2 


- - kp + - q - E 2 

2mab V 


1/2 


(7) 


1/2 


I [ rUab . „ 

-- kp +- q - E 2 

2mab V "ia 


1/2 


( 8 ) 


for E^f,^ < E 2 < Ef^ and hnally 

1 

E 2 


2m, 


ab 


mab 


m„ 


1/2 


Eo 


2m, 


ab 


, , mab 

kp H- q 

mn. 


1/2 


(9) 


for E 2 > E+. 

A comparison between the absolute value of the matrix elements with the analytic 
extension in (|6]), with E 2 ^ E 2 + ie in ([2]) and numerically calculated from ([1]) is shown 
in hgure O The analytic extension becomes more relevant as E 2 approaches E'^j^. For 
E 2 < Eff^ the analytic form in ([2]) nicely reproduces the numerical calculation. 

The correction introduced by the analytic extension can be seen in the spectral 


molecular function, dehned as Amoi{E,q) = —2A T{E,q). Figures 3(a) and 3(c) show 


the spectral molecular function respectively for = 0.1 and 2.5 calculated with (|2]), 
where a small imaginary part is added to the two-body energy. On the other hand, 
hgures 3(b) and 3(d) are calculated with ([6]), where the analytic extension of the matrix 
elements has been properly taken into account. Notice that the use of the proper 


expression slightly changes the contrast in the figures. Figures 3(a) and 3(c) have been 

obtained in |3l] and in order to straightforwardly compare £gure[3]to the previous work, 

calculations here were made replacing the binding energy in (|2]) and (jO]) by the total 

two-body energy measured in terms of the Fermi energy Ep plus the chemical potential 

2 

of the impurity (/ia), specihcally E 2 ^ E 2 — those equations. 
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Figure 2. Comparison between the absolute value of the matrix elements with the 
analytic extension in (|6]), with E 2 ^ E 2 + ie in ([2]) and numerically calculated from 
(P) as function of E 2 /\Eab\ for nia = mb, q/\/ma\Eab\ = 2 and kp/\Jma\Eab\ = 3 . 
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(a) =0.1, without analytic extension. 



0 0.5 1 1.5 2 2.5 3 


q/kp 

(c) = 2.5, without analytic extension. 

Figure 3. Spectral molecular function Amoi{E,q) for = 0.1 (top) and = 2.5 
(bottom). The left-hand-side is the result of calculations with ([2]) and is the same 
as presented in [34]. The right-hand-side is calculated with ([6]), where the analytic 
extension of the matrix elements has been properly taken into account. Notice how 
the contrast change from left to right and that the corrected expression for the matrix 
elements is more relevant in the large binding limit. 
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A way of taking into account fluctuations in the background medium is by introducing 
the self-energy of the impurity, which considers particle-hole excitations on top of the 
Fermi sea. The full propagator for the self-energy of the impurity relates to the bare 
one through the Dyson equation 


G^\E,p) = E 


p 


2m„ 


K(E,p). 


( 10 ) 


In the self-consistent way, the transition operator and the self-energy obey a set of 
non-linear coupled equations. This problem is handled by solving this set of non-linear 
coupled equations iteratively, which is equivalent to a Ward-Luttinger approach [41]. 
The zero-order transition operator is the bare one, given in ([S]) and the zero-order self¬ 
energy is a function of the bare transition operator. The first-order transition operator 
is a function of the zero-order self-energy and so on 
This procedure leads to the set of coupled equations 


see 


Appendix A for the details). 


e:(e,,a) = 


iTpi T„ IS, + S 


^n+1 


(-^2, q) — ^ 
f d'^k- 


Pb<kF 

-1 


pI 


’ 2mfe 


+ Ea ) , n > 0 


e{k- kp) 


-Eo 


2m„, 


+ S - b - ’ 


(11) 


( 12 ) 


where Tq^{E 2 , q) is given in (|6]), iria is the mas of the impurity and ruh the mass of the 
fermion. 

The problem of one single impurity on top of a Fermi sea has been successfully 
solved with non-self-consistent methods (see, e.g. [Ml Ea ESI S2l 03]). However, we 
found the self-consistent one more suitable to handle in the investigation of the three- 
body problem of two impurities immersed in a Fermi sea, whose results are presented 
in the following sections. This choice is supported since the self-consistent method 
describes well the two-body system when particle-hole fluctuations are taken into 
account. For instance, calculations done with the zero-order expression of the self¬ 
energy ([TTl with n = 0) for rua = mb reproduce previous results in the literature. Using 
(na and ([T^ we found that the energy of the attractive and repulsive branches of the 
polaron, which are solutions from E = / (flma) + 3? [U°(U, q)], as well as the weights 

of the branches (Z = {1 — dp {Ti^{E, q))]} ) and the polaronic spectral function 

{Aa{E,q) = —23[U -\-ie — are the same as found in [391134] . 

In [361112] the molecular regime, namely the region where the energy of the molecule 
under the effect of particle-hole fluctuations from the medium (dressed molecule) is lower 
than the polaron one, calculated with non-self-consistent methods happens from large 
and negative r] up respectively to -0.97 or -0.95, above which the polaron energy will 
always be below the molecular energy. We define rj = \n{kpa 2 D) = 

020 is the two-dimensional scattering length which fulfills Eab = iflmahoiD) ^ 

Ep = kf/2mb is the Fermi energy. This transition between regimes is often called 
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polaron-molecule transition and we fonnd a first-order (n = 0 in fllip and 0121 ) ) polaron- 
molecnle transition to occnr aronnd rj = —0.54, which is below the resnlts obtained 
with variational |36] and diagrammatic Monte Carlo 021 calcnlations, bnt is close to 
the r] ^ —0.60 fonnd in the 2D limit of qnasi-two-dimensional highly polarized Fermi 
gases [43] . 

The different behavior fonnd between onr method and previons stndies in the small 
window —0.96(±0.01) < rj < —0.54 [361112] or —0.60 < rj < —0.54 [43] is not important 
for onr pnrpose of nsing the self-consistent approach in the three-body calcnlations. 
The actnal point of crossing is not the relevant information since experiments [31] have 
shown that the crossing of the energies of the polaron and the molecnle as fnnction of r] 
happens at a very shallow angle and that the two-body energy spectrnm is not affected 
by the polaron-molecnle transition (see [31] for details). 

In view of that, it is more important that the model reprodnces this transition 
than the actnal point of occnrrence. Solving (ITT]) and (14^ with n = 0 we find that the 
energy spectrnm of the polaron and the dressed molecnle are very similar to each other 
and that the crossing between the energies does happen at a shallow angle, as shown in 
fignrelH Here a difference between the self-consistent method and the non-self-consistent 
one can be seem. Comparing fignre |4] to [43] , we note that the polaronic spectrnm is 
exactly the same, while the spectrnm of the dressed molecnle has a small qnantitative 
difference. However, as explained above, since onr model captnres the correct behavior 
of the two-body system, this small qnantitative difference is not a problem. As an 
extra argnment for nsing the self-consistent method in our analysis, notice that our 
three-body results when particle-hole fluctuations are taken into account lie mostly in 
1 ] > —0.54 (see hgures [HI and fTOl) . from where the energy of the polaron is always 
lower than the molecular one and it is the same as calculated with non-self-consistent 


methods 

It was found in m, with diagrammatic Monte Carlo calculation, that corrections 
beyond one particle-hole in the medium does not contribute significantly to the energy 
spectrum of the polaron and the dressed molecule. In view of that, the convergence of 
the iterative solution of ffTT]) and ffT 2 l) is studied. Yl^{Ea,0) and t~^i{E 2 , 0 ) are shown 
as function of E 2 /Eab for different 77 ’s, where the strong binding-limit corresponds large 
negative rj and the weak binding-limit to large positive rj. 

Notice that the self- 


Results for large binding {rj = —2) are shown in figure 5(a) 


energy converges very fast as E^(ii^ 2 ) and IIKE 2 ) are almost indistinguishable from each 
other and both of them differ just slightly from I1^{E2). As expected, the self-energy 
correction is almost negligible in this limit, which is confirmed by noticing that Tq^{E 2 ), 
Tf^{E 2 ) and Tf^{E 2 ) are practically identical. 

Increasing kp and moving towards the small binding limit, a few iterations are still 
needed for the self-energy to converge. For 77 = 0, for instance, the difference from 
E^(i? 2 ) to S^(F^ 2 ) is clearly much smaller than the difference from S^(F^ 2 ) to T,^{E 2 ), as 
shown in figure 5(b) The difference is that now the self-energy influences the two-body 
system, rendering it more bound. There is a big difference between the pole position 
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Figure 4. The energy spectrum of the polaron (attractive branch) and the dressed 
molecule as function of the parameter rj. Notice that the crossing between lines happens 
at a shallow angle and that both spectra are similar in the region —1.0 < y < —0.54, 
from where the polaron has always a lower energy. 


of Tq ^{E 2 ) and Ti but this difference is already much smaller from ^(-E' 2 ) to 

Tp(E,). 

The scenario keeps this pattern for rj > 0 and hgures showing details are not 
presented here. The results in hgure |5] provide our argument for working with Tf^(E 2 ) 
and S°(i? 2 ) in the three-body calculations below. The computational cost of going 
beyond this are considerable and we do not expect neither qualitative nor large 
quantitative changes. 

One great simplihcation comes by noticing that the function S is almost 
independent of the angles in its argument, namely 



d9 -^^ 

ME2, q,k) + ^^ cose + S-(/2(E2, q, k), q- k) 



de - 

/i(-^2, q, k) -\- 


1 

^cos0 + E-(/2(E2,g,fc),g2 + F) ’ 


(13) 


where ME2, q, k) = -E2 + - te and /2(^2, q, k) = E2 + 2(j!+m,) “ &,■ 

The agreement between the both sides of flT3l) is in general better than 0.1%. The 
approximation in flT^ is used to calculate and E 2 from flTT|) and ([12]), as shown in 
hgure [5l 
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(a) 77 = -2 (b) 77 = 0 

Figure 5. The transition operator and the self-energy EJJ as functions of E 2 /Eab 
for n = 0 , 1,2 and q/kp = 0 . 


4. Thresholds of the three-body system 


We now treat the problem of two identical impurities of mass roa immersed in a Fermi 
sea background and interacting with a fermion of mass rub that belongs to the sea. 

We assume that the total angular momentum is zero for both interacting and non¬ 
interacting impurities. As for the two-body system, the interaction between particles 
is assumed to be attractive s—wave zero-range potentials and now the center-of-mass 
(CM) of the three-body systems is considered at rest in the frame of the Fermi sea. 
The coupled homogeneous integral equations for the bound state of an aab system are 
derived in detail in [19]. Here we have a Fermi sea and need to take Pauli blocking 
into account. We do so by inserting appropriate ©-functions in the coupled integral 
equations. The Faddeev components describing the system of two impurities immersed 
in the Fermi sea are given by 


fb (<h) = E3 




aa,b 


fa (^) = Tab E3 - 


2m, 

2.mab,a 

d‘^k — 


dfk 


0 (g - kp) fa (k 






d^k- 


2mab 2maa rn, 

e{k-kF)fb(k 


—fc ■ q 


E^ 




3 2maa 2mab rua 

Q {\q + k\- kp) fa (k 


—■ q 


— + -^k-q 




(14) 


(15) 


where the three-body reduced masses are maa,b = 2,mamb/{2ma -\- mb) and mab,a = 
nT'ai.fna + -b mb). The ab and aa transition amplitudes are respectively given 

in (I2|) and (I3|) and are calculated for energies of the corresponding subsystems within 
the aab system. 

Before seeing how the presence of the Fermi sea affects the three-body states, it is 
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important to understand how the two- and three-body thresholds move with the Fermi 
momentum. As the a particles are not directly affected by the Fermi sea, the transition 
operator of this subsystem is given in ([3]). The two-body breakup due to the interacting 
identical a particles where aab —>■ aa + b, is defined as the minimum energy which 

satishes {E^Jf — (f /2maa,b) = 0. The assumption that the three-body system is at 
rest in the Fermi sea frame implies that the total momentum of the pair is restricted 
by the momentum of the fermion (6 particle), i.e., q > kp. The momentum which 
minimizes the energy is kp and Eff moves with kp as 


-rpth _ 

Eaa - 


\Eaa \ + 


kf 


(16) 


In 2D we set Eaa = 0 when the a particles are not interacting with each other 
In this case the two-body transition operator in flTT)) diverges and decouples the 
homogeneous integral equations o and UHl. Besides, the aa + h decay channel 
is not available and Eff need not be considered in the three-body calculations. 

In the same way, the two-body breakup due to the interaction between the ab 
particles {E^^), where aab ab a, is dehned as the minimum energy which satishes 

i^ab ~ ^/‘^'^ab,a) = 0. As the b particle does interact with the Fermi sea, the 
transition operator is given in ([2]) and the expression which satishes {E 2 ) = 0 is 


presented in (jl]). Replacing E 2 E\ 


th 




in this equation allows us to calculate the 


minimum three-body energy that gives a pole in the two-body T-matrix. Notice that 
there is a subtlety in this part since in cases where kp > ^y8mab\Eab\, the minimum 
energy which dehnes the two-body breakup does not correspond to an ab bound state 
(see hgure 1(b)). When the two-body subsystem is virtual (unbound), the three-body 
system will only be stable while its energy is below the lowest two-body threshold. 
This condition is found by calculating the expression which minimizes the energy in the 
expression for Eff^ when Ef^^ —)■ E^Jf — ^ . Collecting everything, the dependence of 

E^Jf with kp reads 



— Eab 

+ 

kl 

2mab 

— \Eab\ 


^l\Eab\ 


^^ab^a^^ab 

kl 1 


mab\ 


2maJ 


>2ml\E, 


ab 


ki 


^ab,amlb 


kp < k* (17a) 

k* <kp< k* (^1 + 

kp >k*(l + (17c) 


where k* = {2ml\Eab\/mab,af^‘^. 

The result in (fTTD is illustrated in hgure El The three-body system has to be below 
either the ab subsystem threshold or the two-body continuum when the hrst one is not 
available. 

Finally, the three-body breakup (Ef^), where aab a + a + b, is dehned as the 
minimum energy which renders one of the denominators on the right-hand-side of (ITT)) 
and fllSp zero. For kp = 0 this energy is E^^^ = 0 and a hnite kp increases the breakup 
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Figure 6. Two-body total energy, E 2 /\Eab\, as function of q/ymif\Efb\. For 
kp < ^y8mab\Eab\ the relevant threshold for the three-body aab system is the bound 
ab subsystem, while for kp > y^8mab\Eab\ it is the unbound subsystem. 


value. Searching for the minimum energy which makes the denominator vanish on the 
hrst term of both o and (na, we have that one of the momenta ^ or fc is allowed 
to vary in the range [0, oo], while the ©—function restrict the other one to [fcir, oo]. In 
both cases, the minimum energy that makes the denominator zero is = kf/2mab- 
In the missing term, the second one on the right-hand-side of flThll . both momenta are 
allowed in the range [0, cxo], but the ©—function restricts this whole term to contribute 
only when q + k>kp (see flB.Tap h The dependence of on kp is therefore 


pth _ 

“ 4m„ 


(18) 


which is always below k‘f/2mab- Equation fflSj) thus defined the three-body threshold 
for the breakup aab a a b. 


5. Interacting impurities 

Each impurity particle, a, interacts with a fermion from the sea by a zero-range potential 
characterized by a two-body binding energy Eat = {2mabalp,) \ forming three-body 
bound states aab. These states can exist even when the impurities are not interacting 
with each other. In both the interacting {Eaa 7 ^ 0) and non-interacting {Eaa = 0) 
cases, the three-body energies increase with increasing Fermi momentum kp. Since the 
two- and three-body thresholds are also moving up as kp increases, three-body bound 
states are found below the Fermi energy even with positive energies. As an example, we 
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show in figure [7] the case where the two impurities have the same mass as the fermion 
= mb). 

The two impurities can, in principle, interact with each other with any energy. We 
set Eaa = Eab, which for identical masses and kp = 0 leads to the well-known case of 
just one two-body state and two three-body bound states with energies E^'^ = 16.52ii^2 
and = 1.27T^2 [Ml US]. As seen in figure [71 the energies and thresholds increase 
with increasing kp. For Eaa = Eab, the decay channel aab a -\- ab is always closed. 
The three-body threshold is highest for kp = 0 and moves up at the smallest rate. It 
goes below Eff when kp/^Jma\Eab\ > \/2. It is interesting to note that the first excited 
state decays into an atom and a dimer when kp/ma\Eab\ ~ 1-1, while the ground 
state is breaking up into three free atoms for kp/^Jma\Eab\ ~ 5.1. 

In view of what was discussed in the previous section, for kp > k* = 
{2mf\Eab\/mab,a)^^‘^ three-body states of the two interacting impurities, which in 2D 
forms a L=0 bound state, and a fermion from the sea are supported even when two of 
the subsystems are in the virtual state, as shown in figure [71 



Figure 7. Three-body energy, E^/\Eab\, as function of the Fermi momentum, 
kp/y^ma\Eab\, for nia = mb and Eaa = Eab The f marks the point kp = k*. For 
kp > k* the two a6-subsystems are unbound. 

The results in figure [71 indicate that cold atomic experiments may be able to deepen 
the knowledge of 2D three-body systems by looking at two impurities in the presence 
of a Fermi sea with which the impurities interact. The Fermi momentum is given by 
density of the sea, which is a controllable parameter. Then, it is possible to change it 
and measure the number of atoms lost in the trap, which must have a peak when the 
three-body system is close to either the two- or three-body continuum |6]. Different 
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bound states may decay into different systems, as seen in figure [7] and a larger number 
of bound states can be reached in highly asymmetric mass systems {mb/rna <C 1) [22] , 
In order to give some insight for experiments, it is possible to systematically extend 
the procedure used in the investigation of the symmetric mass case and study how the 
three-body states vanish for a large range of mass-ratios. The final result is summarized 
in a mass versus kp diagram, shown in figure [8] for Eaa = Eah- The Arabic numerals 
indicate the number of bound states in each region, where 0 means that no bound state 
was found, 1 indicates that only the ground state is available and so on. The central 
black-line divides the plot in two main regions where the three-body states go into 
either a dimer plus atom or three-atom continuum. Taking as example two ^^^Cs atoms 
immersed in a sea of ®Li fermions, the mass-ratio is mb/rua = 0.045 and four bound 
states are present for both Eaa = 0 and Eaa = Eab when kp = 0 [Ml El]. It is shown in 
figure [HI that the four states for Eaa = Eab disappears at r; 1.96 (ground), rj ^ 0.59 
(first), 7] ^ 0.058 (second) and rj ^ —1.42 (third), respectively. Furthermore, notice 
that the three deeper states decay into three atoms while the highest one decays into 
atom plus dimer. 



Figure 8 . Bound state diagram for Eaa = Eab as function of the interaction strength 
rj. The Arabic numerals indicate the number of bound states in each region. The “-I-” 
sign in 4-|- indicates that more than four bound states can be found in that region. 
The central black-line divides the plot in two main regions where the decay channel is 
either aab b aa or aab a + a b. 

The result shown in figure[8]is, to the best of our knowledge, the first one to consider 
the problem of two impurities immersed in a Fermi sea in 2D when the impurities 
are allowed to interact with each other. These results were achieved without taking 
into account particle-hole fluctuations in the Fermi sea. We now consider fluctuation 
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corrections starting with two non-interacting impurities. 

6. Self-energy corrections to the three-body system 

Here we study the three-body system for two non-interacting impurities and then we 
see how the fluctuations in the Fermi sea affect the results. Taking Eaa = 0 leads to 
fbio) = 0 in flTT|) and fflSj) . simplifying the calculation. Furthermore, for = mb and 
kp = 0 only one bound state is supported |16] as shown in flgure|9l Without considering 
the particle-hole fluctuations in the Fermi sea, the behavior of the three-body energy 
and of the thresholds is similar to the case where the impurities are allowed to interact. 
For kp/> 0, the three-body threshold, which is the highest for kp = 0, 
increases at the smallest rate and for kp/ ma \ Eab \ > 1.155 it becomes lower than the 
ab threshold, meaning that the decay channel aab a ab is not available any more 
for the three-body system. When kp/ma\Eah\ ~ 1.7 the three-body state touches the 
continuum and the system decays into three free atoms. 

As in the case of interacting impurities, for kp > k* = {2ml\Eab\/mab,a) ' the 
three-body ground state is supported even when the two afe-subsystems are in a virtual 
state. Since the last pair is also unbound (non-interacting impurities), the three-body 
system is bound when the three subsystems are unbound. The behavior of the ground 
state is similar to what has been seen in so-called Borromean systems in 2D [TTl H8] . 
In the present case we have the background medium in the form of a Fermi sea which 
is what makes this behavior possible. One may therefore consider this an example of a 
medium-induced Borromean state in 2D. 

Next we introduce the self-energy correction in the three-body equations. Since 
the two impurities are not interacting, only the second term on the right-hand-side 
of (1151) survives. It is important to emphasize at this point that the bosons in our 
problem are not from a Bose gas, but isolated impurities, or we can say that the average 
distance with a third boson from the gas is very large compared to any size scales in 
the problem (vanishing limit of the density). In this idealized situation, the connected 
three-body T-matrix can in principle contribute only to the two-boson T-matrix with a 
term proportional to the fermion density. If in addition, we take into account that the 
Bose gas density is negligible according to our statement of the physical situation, the 
contribution of the connected three-body T-matrix to the self-energy will be negligible 
and only the non-connected term of the T-matrix with the fermion and boson T-matrix 
will be relevant. Indeed, this term was already taken into account in our evaluation 
of the self-energy of the impurity. In a more general situation, e.g., in the presence 
of a Bose gas, other contributions to the impurity self-energy should be included, and 
the feedback from the full three-body T-matrix should be relevant as the boson density 
increases. In this case, one has to include corrections to the self-energy beyond the 
contribution of the disconnected part of the three-body transition matrix. 

That said, the first correction comes in the two-body sector where the transition 
operator in ([2]) is replaced by the one in (IT^ with n = 0. As we discussed before. 
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the convergence of the two-body T-matrix with n is very fast and the most signihcant 
correction is seen in £gure|5]to be from to justifying the choice of n = 0 in flT^ . 
In order words, the real gain of going to higher n is not worth the extra time spent in 
the calculation. Then, considering the dressed propagator of the impurity, the integral 
equation is written as 

^ \ f dE2 _ d^ _ 

J J 2m E 2 -£-^- ^,(^ 2 , k) + te 


/(^ = -Ti Es - 


X 


0 



E^ — E 2 — Ea 


(g+fc)^ 

2mh 


+ ie ’ 


(19) 


with given in ffT^ . The numerical solution of flT^ gives the corrected three-body 
energy and we need to understand again how the thresholds and energies depends on 
the Fermi momentum. However, since the self-energy correction is being considered, 
both the two- and three-body corrected thresholds are calculated numerically. 

As shown in hgure [5l the inclusion of the self-energy in the propagator of the 
impurity renders the molecule more bound, which lowers the ab threshold in the three- 
body calculation (see hgure [9]). The same happens to the three-body threshold, which 
starts from zero when kp = 0 and goes negative for kp > 0. The three-body state is 
also more bound and its energy as function of kp increases slower than in the previous 
cases, where corrections due to huctuations of the medium were not considered. The 
three-body energy and the thresholds with and without the self-energy correction are 
compared in hgure [H The hnal result is that, although rendering both the two- and 
three-body bound states more bound, the inclusion of the self-energy makes the states 
disappear at a slightly smaller kp, since the changing rates of the thresholds are strongly 
ahected. 

The extension of the result in hgure [9] for another mass ratios gives the diagram 
shown in hgure [101 where the Arabic numeral indicate the number of bound states in 
each region. The continuous lines are from calculations without huctuations in the Fermi 
sea and the discrete points show how the lines move when the self-energy is considered. 
Notice that the correction is always small within the range of the Fermi momentum 
and masses considered, which covers typically experimentally accessible cold atomic gas 
systems. Such correction is expected to be smaller in the case of interacting impurities, 
where the extra attraction would render the ehects of the Fermi sea less relevant. 








mb/nia Es/lEabl 
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Figure 9. Three-body energy, E^/lEabl, as function of the Fermi momentum, 
kp/\/fna\Eab\, for nia = mb and Eaa = 0. Results with and without the self-energy 
correction are compared. 



Figure 10. Bound state diagram for Eaa = 0 as function of the interaction strength 
rj. The Arabic numerals indicate the number of bound states in each region. The “-I-” 
sign in 3-|- indicates that more than three bounds states can be found in that region. 
The central black-line divides the plot in two main regions where the decay channel is 
either aab a + ab or aab a + a + b. 
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We have considered the 2D problem of two identical atomic imparities (either bosonic or 
with distinct internal states), immersed in a backgronnd Fermi sea and interacting with 
a fermion on top of it in 2D. The interactions were modeled by attractive zero-range 
potentials and the Faddeev decomposition was used to write the homogeneous coupled 
integral equations for the three-body bound state. 

The problem of just one impurity propagating in the Fermi sea was hrst considered 
and we have reviewed some well-known results in the literature |32l |33l EH ESI 1^ - 
specihcally the two-body T-matrix in the medium ([2]) and the energy of the molecule 
as function of its total momentum (|1])- The binding energy of the molecular state was 
studied previously as a function of its mass ratio, total momentum, and Fermi energy. 
For some values of these parameters the pole of the T-matrix, which represents a bound 
state, enters through the cut of the two-body continuum and increasing the momentum 
of the molecule it returns from the cut and become a bound state again ISlllQ]. We 
have shown here that this state appears in the second energy sheet and becomes a 
virtual state, as shown in hgure 1(b) Replacing E 2 ^ E 2 + ie in order to consider 


the two-body T-matrix in the medium ([2]) in the scattering region leads to inconsistent 
results in some cases (see hgure [2]). We have shown that the analytic extension of these 
matrix elements, as presented in ([6]), gives the right solution in this case. 

Particle-hole huctuations were then considered and the self-energy of the impurity 
propagating in the sea was self-consistently calculated. The non-linear couped equations 
for the transition operator and the self-energy ffTTj) and ffT^ were solved iteratively and 
we showed that the convergence of both quantities with the number of iterations was 
very fast. Although this method has some disadvantages, we found it suitable to use 
mainly in the three-body calculation. This choice is supported since the self-consistent 
method employed in this work correctly describes the two-body system when particle- 
hole fluctuations are taken into account. We found that the energy of the attractive and 
repulsive branches of the polaron, as well as their weights and the polaronic spectral 
function are the same as found in [39l [SH |13]. The polaron energy is the relevant 
information for rj > —0.54, where most of the effects of the medium on the three-body 
system are calculated (see hgures E] and [lO]) . 

The formalism used here allowed us to study the complex problem of two interacting 
impurities propagating in the Fermi sea. Interestingly, three-body states of the two 
interacting impurities and a fermion of the sea are supported even when two of the 
subsystems are in a virtual state. Importantly, the fate of the three-body state depends 
strongly on the two-body systems. This feature can be used to identify the three-body 
states of two impurities immersed in a Fermi sea in measurements of cold atoms lost 
from a trap. 

The complexity of the integral equations when the impurities are allowed to interact 
makes the inclusion of the self-energy correction very hard to be implemented. However, 
if the impurities are not interacting, the integral equations simplify and the effect of the 
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correction can be studied. First of all, the three-body bound states are still allowed 
when the impurity-fermion subsystems are in a virtual state. Since we have two non¬ 
interacting impurities, the three-body systems are bound when the three two-body 
subsystems are unbound, which can be interpreted as a medium-induced Borromean 
state in 2D. The fate of the three-body states once more strongly depends on the two- 
body subsystems. The inclusion of the self-energy of the impurity drastically affects the 
behavior of the subsystems and of the three-body states, as it can be seen in figure [9l 
but the superposition of all the effects together leads to just a small correction in the 
hnal results, as shown in figure [101 We expect that the correction would be smaller in 
the case of interacting impurities, since the effects of the Fermi sea would decrease as 
there is more attraction in the system. 

Three-body bound states were also found when a single impurity is immersed in 
a Fermi sea [36] • The different particle can be in a polaronic state or bind one or two 
fermions from the sea. It is interesting to understand what scenario is more favorable 
when another impurity is brought into the game. Would each impurity bind to two 
fermions or the two impurities bind to one fermion? Valuable information can also be 
gained by considering two polarons interacting with each other as a four-body system 
of two fermions and two impurities, using techniques similar to the ones employed 
in |l9]. Another interesting direction for future investigations would be to study how 
the presence of the Fermi sea would affect the momentum distribution of the three- 
body systems by calculating the measurable two- and three-body contacts through 
an extension of the techniques described in [50]. Another interesting point is the 
understanding of how the signature of the Efimov effect in 3D three-body systems under 
the influence of a Fermi sea [51] would disappear and connect to the result presented 
here through a change in dimensionality, similarly to what was done for three-identical 
bosons in [52] [53] . 
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Appendix A. Equation for the self-energy 

The bare and full propagators and the self-energy of the impurity of mass ma immersed 

in a sea of fermions relates to each other through the Dyson equation Gf}{E^p) = 
2 

E — - T,a{E,p). In the self-consistent way, the transition operator is written as 

f{Ea,Pa,E'^,p'^]Eb,Pb,Elpffj =iX + {iXfl{E,q) + {iXfl^{E,q) + ... 

(A.l) 


A-i-/(E,g) ’ 
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where the total energy and total momentum are respectively E = Ea-\-Ei, and q = Pa+Pb- 
The integral I{E,q) reads 


I{E,q) = 


27r 


dfk- 


Q{k- kp) 


E-E' 


{q-ky 

2m„ 


Ta{E - El,q 


fc2 


k) + ie Efj 2 mb 


and the expression for the self-energy is found to be 

0 {kp - Pb) 


T.a{Ea,Pa) = - d^Pb 


^ ^ ~ HEa + ^^,Pa -\-Pb) 

The expressions for the transition amplitude and the self-energy, given respectively 
in flA.ip and flA.3D . form a system of coupled equations. We solve this system iteratively 
by hrstly setting the zeroth-order of the self-energy in flA.Sp as 


(A. 2 ) 


(A.3) 


KiEa,Pa) = 


>Pb<kF 


d^Pb To (\Pb+Pa\, 


pI 

2mi 


+ E, 


(A.4) 


where the zeroth-order transition operator tq is the one given in ([6]) and is independent 
of Sq. Then, is introduced back in flA.ip . leading to the hrst-order operator ri, which 
requires (1A.2I1 to be solved. Performing the integral on the energy in (lA.2p is not trivial 
and an analysis of the poles has to be done. In the complex plane of E'^^ the poles of 
I{E, q) are 

(A.5) 


El = 


2mb 


— le 




(A. 6 ) 


The pole in flA.Sp is always on the lower half complex plane of Ej^, as the only imaginary 
part comes from the small contribution ie and, for = 0, the pole in flA.bp is on the 
upper half plane. The self-energy is complex and contributes to the location of the pole. 
The study of its behavior and the results for the polaron spectral function in [3l] show 
that Q'(Sa) < 0 for any {k, q, E, E'f), ensuring that the pole in flA.bp is on the upper-half 
complex plane of E^ even for 7 ^ 0 . 

It is also necessary to study the analyticity of r and E^ as function of E and E^. 
The coupled equations flA.2l) and flA.3p indicates that both r and E^ must be analytic in 
the same region. From ([1]) is possible to see that r, and consequently E, are analytic on 
the upper half complex plane of E and therefore on the lower half complex plane of Eb- 
As the pole in flA.Sp is in the semi-plane where the functions r and E^ are analytic, the 
contour is closed through the lower half complex plane of E'^ (counterclockwise). The 
integral in flA.2ll is calculated using the residues theorem and leads to the transition 
operator (see flA.ip i 

Tf^{E2,q) = 

f d^k Q{k- kp) 


-E 2 -l- 
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(A.7) 


le 










































Three-body bound states of two bosonic impurities 


23 


where A ^ is related to the energy of the pair in the vacuum (see ([T])). The set of 
decoupled equations for the self-energy and the transition operator are 


S”(^a,Pa) = / d‘^Pb Tn (\pb + Pa\ , 7?^ 

Jpb<kF V 

T'n+\(^2,g) = A-^ 


+ Ea] , n >0 


- / d^k- 


Q{k- kp) 


J -E 2 + ^ + + T2. COS0 + S-(g - k, E^ + - ie 

^ 2mab rria ' aV'i ) z i 2(ma+mb) 2mb' 

where tq is given in ([6]) and is independent of Sq. 


(A,8) 


- .(A.9) 


Appendix B. Angular part of the integral Equation 

Since we are interested in states with total angular momentum zero the spectator 
functions in flIT|) and flT^ are independent of the angle between the momenta q and k. 
The angular dependence of the first term on the right-hand-side of these equations comes 
only from the denominator, as both the spectator {f{q) = /(<?)) and the ©—functions 
are angle-independent. 

However, the ©—function in the second term of flT^ does depend on the angle, 
which makes the integration of such term not as simple as the previous two. The trick 
to solve this integral is to transfer the information of the ©—function to the limits of 
the angular integral. The angular dependence in the ©—function on the second term of 

(USD is 

-|- 2kqcos9 > kp , (B-1) 

which dehnes the angle 6i as the lowest limit of 6 in flB.ip as 


a 

costal = - 




q 


2kq 


(B.2) 


with 0 < 9i < Ti. 

Changing integration limit from 27r to 9i, the angular integral in the second term 
of (irSP reads 


rdi 


d9 


4 tan 


-1 


lath'd-) 


1 -|- a cos 9 


\/l — 


(B.3) 


where the main branch of tan6' has to be considered. Notice that for 9i = ti, flB.3D gives 
the result for the angular integration of the hrst term on the right-hand-side of fflTp and 

(USD- 

The dependence of the ©— function on the angle is eliminated using that | cos^il < 
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2 n d9 Q (\q k\ — kp 
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q k < kp (B.4a) 
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After the angular integration, the homogeneous coupled integral flT4ll and flTdl) read 
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where the function g{q, k, kk) reads 
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